Thin films versus 2D sheets in layered structures: graphene and 2D metallic sheets 
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We study an interface between two media separated by a strictly 2D sheet. We show how the 
amplitude reflection coefficient can be modeled by that for an interface where the 2D sheet has been 
replaced by a film of small but finite thickness. We give the relationship between the 3D dielectric 
function of the thin film and the 2D dielectric function of the sheet. We choose graphene and a 
2D metallic sheet as illustrative examples. This approach turns out to be very useful when treating 
graphene or graphene like sheets in non-planar structures. 
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I. INTRODUCTION 



II. AMPLITUDE REFLECTION COEFFICIENT 
AT AN INTERFACE 



There are no strictly 2D (two-dimensional) systems in 
the real world. However, if the carriers are strongly con- 
fined in one direction they have quantized energy levels 
for one spatial dimension but are free to move in two spa- 
cial dimensions*. Thus the wave vector is a good quan- 
tum number for two dimensions but not for the third. 
Examples are narrow quantum wells and graphene. In 
many situations these systems can be approximated by 
the idealized systems of strictly 2D sheets. Even thicker 
films where the spatial quantization is negligible behave 
as strictly 2D from distances much larger than the film 
thickness. The dielectric functions of a 2D electron gas 2 
and for graphene^— have been derived in the literature. 
The dielectric function for graphene is valid for an ideal- 
ized 2D sheet. If one were to treat graphene as a film of 
finite thickness one could not use the 2D dielectric func- 
tion as it is. It has to be modified into a 3D version. 
The same thing applies to a strictly 2D metal sheet if 
it were to be treated as a film of finite thickness. If a 
metal film of finite thickness with a bulk dielectric func- 
tion were to be treated as a strictly 2D sheet the dielec- 
tric function had to be modified into a 2D version. How 
these modifications can be done is what this paper is all 
about. Examples where these results are useful are for 
optical properties in layered structures and for disper- 
sion interactions (van der Waals and Casimir) in layered 
structures. 

In Sec. [H]we give the amplitude reflection coefficients 
for an interface between two media, in Sec. Ill Al we show 
how these are modified when a 2D sheet is inserted at the 
interface, in Sec. Ill Bl we show how these are modified 
when instead a film of finite thickness is inserted, and 
in Sec. Ill CI we show how a 2D sheet can be modeled 
by a film of finite thickness. In Sees. IIIII and IIVI we 
illustrate this modeling for graphene, and a metal sheet, 
respectively. Before we end with a brief summary and 
conclusion section, Sec. E2 we discuss in Sec. |V] how 
the results for 2D sheets and thin films can be used in 
non-planar structures. 



For the present task we need a geometry consisting of 
two regions and one interface, i\j. For planar structures 
there are two types of mode, transverse magnetic (TM) 
or p-polarized and transverse electric (TE) or s-polarized. 
They have different amplitude reflection coefficients. At 
an interface between medium i and j the TM and TE 
amplitude reflection coefficients for waves impinging from 
the i side are 
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types. If retardation is neglected there is only one 
type and the amplitude reflection coefficient is 
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In the above equations 7$ = y 1 — £j (u) (co/ck) , §i (u>) 
is the dielectric function of medium i, c is the speed of 
light in vacuum, and k is the length of a wave vector 
in the plane of the interface. We have suppressed the 
arguments of the functions in Eqs. ([T][3]). The amplitude 
reflection coefficients and the 7 -functions are functions of 
k and w. The dielectric functions are functions of ui, only, 
i.e., spatial dispersion is neglected. Inclusion of spatial 
dispersion in the bulk dielectric functions is possible 8,9 
but would lead to much higher complexity and negligible 
effects in most cases. 



A. Interface with a strictly 2D sheet 

There are different formulations of electromagnetism 
in the literature. The difference lies in how the conduc- 
tion carriers are treated. In one formulation these car- 
riers are lumped together with the external charges to 
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form the group of free charges. Then only the bound 
charges contribute to the screening. We want to be able 
to treat geometries with metallic regions. Then this for- 
mulation is not suitable. In the formulation that we use 
the conduction carriers are treated on the same footing 
as the bound charges. Thus, both bound and conduction 
charges contribute to the dielectric function. In the two 
formalisms the E and B fields, the true fields, of course 
are the same. However, the auxiliary fields the D and H 
fields are different. To indicate that we use this alterna- 
tive formulation we put a tilde above the D and H fields 
and also above the dielectric functions. See Ref.[6( for a 
fuller discussion on this topic. 

Within this formalism the standard boundary condi- 
tions, used to derive the reflection coefficients, are that 
in absence of external charge and current densities at 
an interface the tangential components of the E and H 
fields and the normal components of the D and B fields 
are all continuous across the interface. The sources to the 
fields in Maxwell's equations are the external charge and 
current densities. In the boundary conditions any dis- 
continuities in the normal component of the D fields and 
tangential component of the H fields are caused by exter- 
nal surface charge densities and external surface current 
densities, respectively. 

The amplitude reflection coefficient gets modified if 
there is a 2D layer at the interface. We treats the 2D 
layer at the interface as external to our system. The 
modified amplitude reflection coefficient for a TM mode 



and the dielectric function is 
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where the polarizability, a", of the 2D sheet is obtained 
from the dynamical conductivity, a", 
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and the dielectric function is 

e 11 (k,cj) = 1 W (k,uj) . 
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For TM modes the tangential component of the electric 
field, which will induce the external current, is parallel to 
k, so the longitudinal 2D dielectric function of the sheet 
enters. The bound charges in the 2D sheet also contribute 
to the dynamical conductivity and the polarizability. 

The modified amplitude reflection coefficient for a TE 
mode is£ 
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where the polarizability, a , of the 2D sheet is obtained 
from the dynamical conductivity, a 1 - , 
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For a TE wave the electric field is perpendicular to k, so 
the transverse 2D dielectric function of the sheet enters. 
The bound charges in the 2D sheet also contribute to the 
dynamical conductivity and the polarizability. 

If retardation is neglected there is only one mode type 
and the amplitude reflection coefficient is 
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Now we have in Eqs. (|4]), ([7]), and (fTO)) the amplitude 
reflection coefficients for an interface between two media 
with a 2D sheet sandwiched in between. To be noted is 
that spatial dispersion of the 2D sheet can be included 
without any complications. This spatial dispersion has 
furthermore important effects as we will show later. In 
next section we will show the corresponding results when 
instead of a 2D sheet we have a thin film sandwiched 
between the two media. 



B. Interface with a film of finite thickness 

For the present task we need a geometry consisting of 
three regions and two interfaces, i\j\k. For this composite 
interface the amplitude reflection coefficient for a wave 
impinging from the i side is n i 12 
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where dj is the thickness of the film j. This expression 
is valid for TM- and TE-modes when retardation is in- 
cluded and also for the modes when retardation is ne- 
glected. The appropriate amplitude reflection coefficient 
from Eqs. (jT][3j) should be used in the expression on the 
right hand side. In next section we show how the effect 
of a strictly 2D sheet can be modeled by a film of finite 
thickness with the proper choice of 3D dielectric function. 



C. Simulating a strictly 2D film with a film of 
finite thickness 

Let us study a 2D-sheet placed in a time varying elec- 
tric field, E, in the plane of the film. There will be a 
surface current density, K = a 2D E. This current den- 
sity has contributions also from bound charges. Since we 
want to treat this 2D sheet as a thin film of finite thick- 
ness, 5, we let this current be spread evenly through the 
thickness of the film. The volume current density, j, is 
then j = K/<5, and since j = er 3J:) E it follows that 
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FIG. 1. (Color online) The attractive nonretardectf interac- 
tion energy per unit area between two graphene sheets as 
function of separation. 



Now, since 
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we find that 
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Thus in problems with 2D sheets one can treat the sheets 
as thin 3D films where the 3D polarizability above is used. 
To check if this is a reasonable approach we insert the 
expressions in Eqs. into Eq. (fTTj) where now dj is S 

and ij = 1 + 2a 2D /k8. If we now let the film thickness, 
S, go towards zero we reproduce the results in Eqs. (j4j), 
([7|), and (fT0|) . Thus in the limit the model is exact. 



III. GRAPHENE AS AN ILLUSTRATIVE 
EXAMPLE 



v. If we now treat the graphene sheet as a thin film of 
thickness S the polarizability of the film material should 
be chosen as 
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and on the imaginary frequency axis it is 
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The van der Waals (vdW) and Casimir interactions can 
be derived in many different ways. One way is to derive 
the interaction in terms of the electromagnetic normal 
modesii of the system. For planar structures the inter- 
action energy per unit area can be written as^i 
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where (wfc) = is the condition for electromagnetic 
normal modes. For the present geometry the mode con- 
dition function is 



fk = 1 - e- 2kd (r 121 ) 



(20) 



where the index 1 stands for vacuum and the index 2 
for the film material. Using Eq. (fTTj) with the proper 
functions for our problem inserted we get 
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For strictly 2D sheets the corresponding mode condition 
function is 



We will now calculate the Casimir interaction energy 
between two free standing undoped graphene sheets in 
vacuum. To make it as simple as possible we neglect 
retardation and perform the calculations for zero tem- 
perature. Retardation effects are actually very small in 
graphene&i£. In a general point, z, in the complex fre- 
quency plane, away from the real axis the polarizability 
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where v is the carrier velocity which is a constant in 
graphene (E = ±hvk), and g represents the degener- 
acy parameter with the value of 4 (a factor of 2 for spin 
and a factor of 2 for the cone degeneracy). In the numer- 
ical calculations we use the valued 8.73723 x 10 5 m/s for 
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The result of Eq. (fT9j) with the mode condition function 
from Eq. (|2"31 is shown in Fig. [T] With the particular 
screening in graphene it turns out that a 2D (k/X, icj/X) — 
a(k,iu>) and the separation dependence of the nonre- 
tarded interaction becomes very simple. A change in 
dummy variables removes the only d in the integrand 
and produces the factor of d~ 3 in front of the integral. 
The result is a straight line in Fig.Q] 

When we treat the graphene sheets as thin films of 
thickness 6 the momentum, k, enters in more places in 
the integrand with the effect that the result has a more 
complicated dependence on d. However in those addi- 
tional places k always enters in the combination kS. This 
means that there will be a universal correction factor, 
from treating the sheet as a film, that depends on d/5. 
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FIG. 2. (Color online) The ratio between the film result and 
the 2D sheet result for the interaction energy between two 
graphene sheets as function of the ratio between the separa- 
tion and the film thickness. 
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FIG. 3. (Color online) The attractive nonretarded interaction 
energy per unit area between two metal films as function of 
separation. They all have the projected 2D carrier density 



This correction factor is shown in Fig. [21 One would 
expect the interaction to decrease for small separations 
since all matter in the two films is further apart than d, 
see the inset in Fig. [5] This reduction is found but then 
there is an unexpected over shoot at larger distances with 
a maximum of 19% at around 86. 



IV. STRICTLY 2D METALLIC SHEETS AS AN 
ILLUSTRATIVE EXAMPLE 



For 2D metallic sheets we use the zero-temperature 
RPA expression: 
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The result for the vdW Casimir interaction energy is 
shown in Fig. [3] for both strictly 2D metal sheets and for 
a series of films with different film thickness. For large 
separations the interaction follows a fractional power 
law-i£~— . At the small separation end of the figure the 
interaction weakens due to spatial dispersion. 

In the metal film case we cannot produce a universal 
curve for the correction factor. Instead we give in Fig. 
131 the correction factor for a series of films with different 
film thickness. However for larger film thickness the re- 
sults approach a universal curve. We see that the overall 
behavior is similar to that in the graphene case but the 
over shooting is here smaller. 
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FIG. 4. (Color online) The ratio between the film result and 
the 2D sheet result for the interaction energy between two 
metal films as function of the ratio between the separation 
and the film thickness. All films have the projected 2D carrier 
density 1 x 10 13 cm -2 . Note that all curves with thickness of 
5A or larger fall more or less on the same curve. 



V. NON PLANAR STRUCTURES 

In this section we discuss how one may proceed in non- 
planar structures, One may, e.g., have a cylinder or a 
sphere coated by a graphene or graphene-like film. The 
spatial dispersion complicates things. In these structures 
the momentum, k, is no longer a good quantum number 
for the normal modes. However, the problems are often 
dominated by long wavelengths. The 3D polarizability 
for a graphene film in the long wavelength limit is 
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Fortunately the wave number is now absent from the ex- 
pression and nothing hinders the use of this expression 
in non-planar structures. 
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For the metal film in the long wavelength limit 



VI. SUMMARY AND CONCLUSIONS 



and hence 
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i.e., the ordinary Drude result for a 3D metal. Once 
again the wave number is absent from the expression and 
nothing hinders the use of this expression in non-planar 
structures. 



We have shown a way to modify the dielectric func- 
tion when a film of finite thickness is used to simulate 
a strictly 2D sheet or vice versa. We have used both a 
graphene sheet and a strictly 2D metal sheet as illustra- 
tions. We further pointed out how to proceed in the case 
of non-planar structures as to avoid problems caused by 
spatial dispersion. 
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